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Abstract. This work is focused on the doubly nonlinear equation 

d u u + d xxxx u + (p - \\d x u\\ 2 L 2 (01) )d xx u + d t u + k 2 u + = f, 

whose solutions represent the bending motion of an extensible, elastic bridge suspended 
by continuously distributed cables which are flexible and elastic with stiffness k 2 . When 
the ends are pinned, longterm dynamics is scrutinized for arbitrary values of axial load 
p and stiffness k 2 . For a general external source / we prove the existence of bounded 
absorbing sets. When / is time- independent, the related semigroup of solutions is shown 
to possess the global attractor of optimal regularity and its characterization is given in 
terms of the steady states of the problem. 



1. Introduction 

1.1. The model equation. In this paper, we scrutinize the longtime behavior of a non- 
linear evolution problem describing the damped oscillations of an extensible elastic bridge 
of unitary natural length suspended by means of flexible and elastic cables. The model 
equation ruling its dynamics can be derived from the standard modeling procedure, which 
relies on the basic assumptions of continuous distribution of the stays' stiffness along the 
girder and of the dominant truss behavior of the bridge (see, for instance, [21]). 

In the pioneer papers by McKenna and coworkers (see [22, 24, 25]), the dynamics of a 
suspension bridge is given by the well known damped equation 

(1.1) d tt u + d xxxx u + d t u + k 2 u + = f, 

where u = u(x,t) : [0, 1] x R — > R, accounts for the downward deflection of the bridge in 
the vertical plane, and u + stands for its positive part, namely, 



;i.2) 



if u > 0, 
if u < 0. 



Our model is derived here by taking into account the midplane stretching of the road bed 
due to its elongation. As a consequence a geometric nonlinearity appears into the bending 
equation. This is achieved by combining the pioneering ideas of Woinowsky-Krieger on 
the extensible elastic beam [30] with equation (1.1). Setting for simplicity all the positive 
structural constants of the bridge equal to 1, we have 

(1.3) d tt u + d xxxx u + (p - \\d x u\\ 2 L 2^ tl ^d xx u + d t u + k 2 u + = f, 
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where / = f(x,t) is the (given) vertical dead load distribution. The term —k 2 u + models 
a restoring force due to the cables, which is different from zero only when they are being 
stretched, and dtu accounts for an external resistant force linearly depending on the 
velocity. The real constant p represents the axial force acting at the ends of the road 
bed of the bridge in the reference configuration. Namely, p is negative when the bridge is 
stretched, positive when compressed. 

As usual, u and d t u are required to satisfy initial conditions as follows 




u(x, 0) = uq(x), iG [0,1], 
d t u(x : 0) = u\(x), iG [0,1]. 



Concerning the boundary conditions, we consider here the case when both ends of the 
bridge are pinned. Namely, for every t G K, we assume 

(1.5) «(0, t) = u(l, t) = d xx u(0, t) = d xx u{\, t) = 0. 

This is the simpler choice. However, other types of boundary conditions with fixed ends 
are consistent with the extensibility assumption as well; for instance, when both ends are 
clamped, or when one end is clamped and the other one is pinned. We address the reader 
to [18] for a more detailed discussion. Assuming (1.5), the domain of the differential 
operator d xxxx acting on L 2 (0, 1) is 

V(d xxxx ) = {we # 4 (0, 1) : w(0) = w(l) = d xx w(0) = d xx w{l) = 0}. 

This operator is strictly positive selfadjoint with compact inverse, and its discrete spec- 
trum is given by A n = n 4 7r 4 , n E N. Thus, Ai = ir 4 is the smallest eigenvalue. Besides, 
the peculiar relation 



Kyxxxx) Uxx 

holds true, with Dirichlet boundary conditions and 

V(-d xx )=H 2 (0,l)nH 1 (0,l). 

Hence, if pinned ends are considered, the initial-boundary value problem (1.3)— (1.5) can 
be described by means of a single operator A = d xxxx , which enters the equation at the 
powers 1 and 1/2. Namely, 

d tt u + Au + d t u - (p - \\u\\l)A 1/2 u + k 2 u + = /, 

where || • ||i is the norm of H^(0, 1). This fact is particularly relevant in the analysis of 
the critical buckling load p c , that is, the magnitude of the compressive axial force p > 
at which buckled stationary states appear. 

As we shall show throughout the paper, this model leads to exact results which are 
rather simple to prove and, however, are capable of capturing the main behavioral dynamic 
characteristics of the bridge. 
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1.2. Earlier contributions. In recent years, an increasing attention was payed to the 
analysis of buckling, vibrations and post-buckling dynamics of nonlinear beam models, 
especially in connection with industrial applications [23, 26] and suspension bridges [2, 3]. 
As far as we know, most of the papers in the literature deal with approximations and 
numerical simulations, and only few works are able to derive exact solutions, at least 
under stationary conditions (see, for instance, [6, 7, 12, 19]). In the sequel, we give a brief 
sketch of earlier contributions on this subject. 

In the fifties, Woinowsky-Krieger [30] proposed to modify the theory of the dynamic 
Euler-Bernoulli beam, assuming a nonlinear dependence of the axial strain on the defor- 
mation gradient. The resulting motion equation, 

(1-6) d tt u + d xxxx u + (p- \\d x u\\ 2 L2(0>1) )d xx u = 0, 

has been considered for hinged ends in the papers [5, 13], with particular reference to 
well-posedness results and to the analysis of the complex structure of equilibria. Adding 
an external viscous damping term d t u to the original conservative model, it becomes 

(1.7) d tt u + d xxxx u + d t u +(p- \\d x u\\ 2 L 2 (0A) )d xx u = 0. 

Stability properties of the unbuckled (trivial) and the buckled stationary states of (1.7) 
have been established in [4, 14] and, more formally, in [27]. In particular, if p < p c , 
the exponential decay of solutions to the trivial equilibrium state has been shown. The 
global dynamics of solutions for a general p has been first tackled in [20] and improved 
in [15], where the existence of a global attractor for (1.7) subject to hinged ends was 
proved relying on the construction of a suitable Lyapunov functional. In [11] previous 
results are extended to a more general form of the nonlinear term by virtue of a suitable 
decomposition of the semigroup first introduced in [18]. 

A different class of problems arises in the study of vibrations of a suspension bridge. The 
dynamic response of suspension bridges is usually analyzed by linearizing the equations 
of motion. When the effects of extensibility of the girder are neglected and the coupling 
with the main cable motion is disregarded, we obtain the well-known Lazer-McKenna 
equation (1.1). Free and forced vibrations in models of this type, both with constant and 
non constant load, have been scrutinized in [3] and [9]. The existence of strong solutions 
and global attractors for (1.1) has been recently obtained in [31]. 

In certain cases Lazer-McKenna's model becomes inadequate and the effects of exten- 
sibility of the girder have to be taken into account. This can be done by introducing into 
the model equation (1.1) a geometric nonlinear term like that appearing in (1.6). Such 
a term is of some importance in the modeling of cable-stayed bridges (see, for instance, 
[21, 29]), where the elastic suspending cables are not vertical and produce a well-defined 
axial compression on the road bed. 

Several studies have been devoted to the nonlinear vibrational analysis of mechanical 
models close to (1.3). Abdel-Ghaffar and Rubin [1, 2] presented a general theory and 
analysis of the nonlinear free coupled vertical-torsional vibrations of suspension bridges. 
They developed approximate solutions by using the method of multiple scales via a per- 
turbation technique. If torsional vibrations are ignored, their model reduces to (1.3). 
Exact solutions to this problem, at least under stationary conditions, have been recently 
exhibited in [19]. 
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1.3. Outline of the paper. In the next Section 2, we formulate an abstract version of 
the problem. We observe that its solutions are generated by a solution operator S(t), 
which turns out to be a strongly continuous semigroup in the autonomous case. The 
existence of an absorbing set for the solution operator S(t) is proved in Section 3 by virtue 
of a Gronwall-type Lemma. Section 4 is focused on the autonomous case and contains 
our main result. Namely, we establish the existence of the regular global attractor for 
a general p. In particular, we prove this by appealing to the existence of a Lyapunov 
functional and without requiring any assumption on the strength of the dissipation term. 
A characterization of the global attractor is given in terms of the steady states of the 
system (1.3)— (1.5). First, we proceed with some preliminary estimates and prove the 
exponential stability of the system provided that the axial force p is smaller than p c . 
Finally, the smoothing property of the semigroup generated by the abstract problem is 
stated via a suitable decomposition first devised in [18]. 



2. The Dynamical System 

In the sequel we recast problem (1.3)-(1.5) into an abstract setting in order to establish 
more general results. 

Let (H, (•, •), || • ||) be a real Hilbert space, and let A : T>(A) d H — > H be a strictly 
positive selfadjoint operator with compact inverse. For r G R, we introduce the scale of 
Hilbert spaces generated by the powers of A 

H r = V(A'/ A ), (u,v) r = (A r /\A r ^v), \\u\\ r = \\A r ^u\\. 

When r = 0, the index r is omitted. The symbol (•, •) will also be used to denote the 
duality product between H r and its dual space H~ r . In particular, we have the compact 
embeddings H r+1 <<= H r , along with the generalized Poincare inequalities 

(2.1) Ai|M|f < \\u\\ 4 r+1 , \fueH r+ \ 

where Ai > is the first eigenvalue of A. Finally, we define the product Hilbert spaces 

U r = H r+2 x H r 

For p G R, we consider the following abstract Cauchy problem on H in the unknown 
variable u = u(t), 



(2.2) 



d tt u + Au + d t u -(p- IMI^A 1 /^ + k 2 u + = fit), t > 0, 
u(0) = M , d t u(0) = Mi . 



Problem (1.3)-(1.5) is just a particular case of the abstract system (2.2), obtained by 
setting H = L 2 (0, 1) and A = d xxxx with the boundary condition (1.5). 
The following well-posedness result holds. 

Proposition 2.1. Assume that f G Ll oc (0,T; H). Then, for all initial data z = (uq,u\) G 
H, problem (2.2) admits a unique solution 

(u(t),d t u(t))eC(0,T;H), 



which continuously depends on the initial data. 
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We omit the proof of this result, which is based on a standard Galerkin approxima- 
tion procedure (see, for istance [4, 5]), together with a slight generalization of the usual 
Gronwall lemma. In particular, the uniform-in-time estimates needed to obtain the global 
existence are exactly the same we use in proving the existence of an absorbing set. 

In light of Proposition 2.1, we define the solution operator 

S(t)eC(U,H), Vt>0, 

as 

z = (u ,«i) i-> S(t)z = (u(t),d t u(t)). 

Besides, for every z G H, the map t h-> S(t)z belongs to C(M + ,"H). Actually, it is a 
standard matter to verify the joint continuity 

(t,z) ^ S(t)z G C(R + x H,H). 

Remark 2.2. In the autonomous case, namely when / is time-independent, the semigroup 
property 

S(t + r) = S(t)S(r) 

holds for all t, r > 0. Thus, S(t) is a strongly continuous semigroup of operators on H 
which continuously depends on the initial data: for any initial data z G T-L, S(t)z is the 
unique weak solution to (2.2), with related norm given by 

CI \ II Il2 II l|2 , II ||2 

t(z) = \\z\\ H = \\u\\ 2 +\\v\\ . 
For any z = (u,v) G H, we define the energy corresponding to z as 

(2.3) E( Z )=S( Z )+ 1 -{\\u\\j-p) 2 + k 2 \\u+\\\ 

and, abusing the notation, we denote E(S(t)z) by E(t) for each given initial data z G %. 
Multiplying the first equation in (2.2) by d t u, because of the relation 

we obtain the energy identity 

(2.4) ^ t E + 2\\d t u\\ 2 = 2(d t u,f). 

In particular, for every T > 0, there exists a positive increasing function Qt such that 

(2.5) E(t) < Qt(E(0)), VtG[0,T]. 

3. The Absorbing Set 

It is well known that the absorbing set gives a first rough estimate of the dissipativity of 
the system. In addition, it is the preliminary step to scrutinize its asymptotic dynamics 
(see, for instance, [28]). Here, due to the joint presence of geometric and cable- response 
nonlinear terms in (2.2), a direct proof of the existence of the absorbing set via explicit 
energy estimates is nontrivial. Indeed, the double nonlinearity cannot be handled by 
means of standard arguments, as either in [24] or in [31]. Dealing with a given time- 
dependent external force / fulfilling suitable translation compactness properties, a direct 
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proof of the existence of an absorbing set is achieved here by means of a generalized 
Gronwall-type lemma devised in [16]. 

An absorbing set for the solution operator S(t) (referred to the initial time t — 0) is a 
bounded set 03% C % with the following property: for every R > 0, there is an entering 
time tR > such that 

|J S(t)z C 

t>t R 

whenever < R. In fact, we are able to establish a more general result. 

Theorem 3.1. Let f e L°°(R + ,f/~) ; and /et d t / fre a translation bounded function in 
L 2 loc (W + ,H- 2 ), that is, 

(3.1) sup / + ||9 t /(r)|| 2 _ 2 dT = M<oo. 

Then, there exists R > with the following property: in correspondence of every R > 0, 
there is t = t (R) > such that 

E(t) < Rq, Vt > t , 

whenever E(0) < R. Both R and t can be explicitly computed. 

We are able to establish Theorem 3.1, leaning on the following Lemma. 

Lemma 3.2 (see Lemma 2.5 in [16]). Let A : R + — > R + be an absolutely continuous 
function satisfying, for some M > 0, e > 0, the differential inequality 

^A(t) + eA(i) <<p(t), 
where tp : 1R + — > 1R + is any locally summable function such that 

rt+l 

sup / v?(r)dr < M. 
t>o Jt 

Then, there exist R\ > and 7 > such that, for every R > 0, it follows that 

A(t) < R u Vt > i? 1/7 (l + 7M)- 1 , 

whenever A(0) < R. Both R x and 7 can be explicitly computed in terms of M ande. 

Proof of Theorem 3.1. Here and in the sequel, we will tacitly use several times the Young 
and the Holder inequalities, besides the usual Sobolev embeddings. The generic positive 
constant C appearing in this proof may depend on p and ||/||l oo (r+,h)- 
On account of (2.4), by means of the functional 

(3.2) C(z) = E(z)-2{u,f), 
we introduce the function 

C(t)=E(t)-2(u(t)J(t)}, 
which satisfies the differential equality 

(3.3) ±c + 2\\d t u\\ 2 = -2(u,dtf). 
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Because of the control 

(3-4) 2\(u,dJ)\< 1 -\\u\\l + 2\\d t f\\ 2 _ 2 , 

we obtain the differential inequality 

^C + 2\\d t u\\ 2 < 1 -E + 2\\d t f\\ 2 _ 2 . 

Next, we consider the auxiliary functional T(z) = (u,v) and, regarding T(S(t)z) as T(t), 
we have 

A T + T + || u ||i + (||«||? - p) 2 +p(\\u\\l -p) + k 2 \\u + \\ 2 - (u, f) = \\d t u\\ 2 . 
Noting that 

^(\\ u \\i - p) 2 + p(\\ u \\l - p) = ^h\\t- 

we are led to 

At + t + 1|« + || u ||i + k 2 \\ u +\\ 2 + \{\\u\\\-p) 2 - (u, f) = \\d tU \\ 2 + \p 2 - 

Precisely, we end up with 

(3.5) l i r + r + lE<l H uf + ±\\ff + l P >. 

Finally, we set 

A(z) = C(z) + T(z) + C, 
where C=^||/|| 2 +2aT + We first observe that A(z) satisfies 

(3.6) \s{z) < -E(z) < A(z) < 2E(z) + c. 

In order to estimate A from below, a straightforward calculation leads to 
A(z) > E(z)-2\(u,f)\-\T(z)\+C 

> E{z) -\\\u\\l- \\\vf -\{\\u\\\- p) 2 - 2\\f\t 2 -±--^= + C 
= \\Hl + \\H 2 + \{\\u\\\ - p) 2 + k 2 \\u+\\ 2 - 2H/II 2 , - ±- - JL + C 

where we take advantage of 

|T(z)| < ||u||||v|| < ^=N|i||u|| < \\M 2 + ^=\\u\\l 

<hv\\ 2 + -^(\\u\\i-p) + bl 



2 V / AT V ' 2j\ x 
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The upper bound for A can be easly achieved as follows 

a(z) < E( Z ) + imi* + inr + l\\u\\t + \\f\\i 2 +^+c 

< 2E{z) + yU\\ 2 + 3^- + C < 2E(z) + c, 

by virtue of 

1 1 11 

(3.7) \T(z)\ < \\u\\\\v\\ < -—lldUM! < \\v\\ 2 + — =||u||? < \\v\\ 2 + -llullf + — — . 

Going back to differential equation and making use of (3.3) and (3.5), the function A(t) = 
A(S(t)z) satisfies the identity 

i A + ^ + l + \\\ u Wl + + ill-Hi + ^> + H 2 = "2<«, dtf) + ^ 

and, as a consequence, we obtain the estimate 

+ £ + i(T + ||u||i + ||^|| 2 + + 4<«, < ^- 

Now, using (3.7) and (3.4), we have 

A A + ^<2||^|| 2 2 + c, 

where c = + y. Thus, by virtue of (3.1) and (3.6), Lemma 3.2 yields 

E(t) < 2A(t) < 2R 1 (M,c). 

□ 

Remark 3.3. If the set of stationary solutions to (2.2) shrinks to a single element, the 
subsequent asymptotic behavior of the system becomes quite simple. Indeed, this occurs 
when p < p c = y/Xi- If this the case, the only trivial solution exists and is exponentially 
stable, as it will be shown in Section 4. The more complex and then attractive situation 
occurs when the set of steady solutions contains a large (possibly infinite) amount of 
elements. To this end, we recall here that the set of the bridge stationary-solutions 
(equilibria) has a very rich structure, even when / = (see [19]). 

4. The Global Attractor 

In the remaining of the paper, we simplify the problem by assuming that the external 
force / is time- independent. In which case, the operator S(t) is a strongly continuous 
semigroup on T-L (see Remark 2.2). Having been proved in Sect. 3 the existence of the 
absorbing set 23, we could then establish here the existence of a global attractor by 
showing that the semigroup S(t) admits a bounded absorbing set in a more regular space 
and that it is uniformly compact for large values of t (see, for instance, [28, Theor. 1.1]). 
In order to obtain asymptotic compactness, the a-contraction method should be employed 
(see [20] for more details). If applied to (2.2), however, such a strategy would need a lot 
of calculations and, what is more, would provide some regularity of the attractor only if 
the dissipation is large enough (see [15], for instance). 
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Noting that in the autonomous case problem (2.2) becomes a gradient system, there is 
a way to overcome these difficulties by using an alternative approach which appeals to the 
existence of a Lyapunov functional in order to prove the existence of a global attractor. 
This technique has been successfully adopted in some recent papers concerning some 
related problems, just as the longterm analysis of the transversal motion of extensible 
viscoelastic [18] and thermoelastic [17] beams. 

We recall that the global attractor A is the unique compact subset of H which is at 
the same time 

(i) attracting: 

lim 6(S(t)<B,A) ->■ 0, 

for every bounded set 23 C "H, where 6 denotes the usual Hausdorff semidistance 
in H; 

(ii) fully invariant: 

S(t)A = A, Wt> 0. 
We address the reader to the books [8, 20, 28] for a detailed presentation of the theory of 
attractors. 

Theorem 4.1. The semigroup S(t) acting on H possesses a connected global attractor 
A bounded in H 2 . Moreover, A coincides with the unstable manifold of the set S of the 
stationary points of S(t), namely, 

A= \z(0) : z is a complete (bounded) trajectory of S(t) : lim \\z(— t) — S\\y, = >. 

Remark 4.2. Due to the regularity and the invariance of A, we observe that S(t)z is a 
strong solution to (2.2) whenever z E A. 

The set S of the bridge equilibria under a vanishing lateral load consists of all the pairs 
(u, 0) e H such that the function u is a weak solution to the equation 

Au-(p- |H|2)A 1/2 m + k 2 u + = 0. 

In particular, u solves the following boundary value problem on the interval [0, 1] 

(d xxxx ii + (bir 2 - \\d x u\\ 2 L2 ( 01) )d xx u + k 2 tt 4 u + = 0, 

\u(0) = u(l) = d xx u(0) = d xx u(l) = 0, 

where we let k = kit 2 , kgR, and p = bir 2 , b G R. It is then apparent that S is bounded 
in H 2 (0, 1) n ^(0, 1) for every 6,kgR. 

When k = 0, a general result has been established in [12] for a class of non-vanishing 
sources. In [6, 7], the same strategy with minor modifications has been applied to problems 
close to (4.1), where the term u + is replaced by u (unyielding ties). 

The set of buckled solutions to problem (4.1) is built up and scrutinized in [19]. In order 
to have a finite number of solutions, we need all the bifurcation values to be distinct. This 
occurrence trivially holds when k — 0, because of the spectral properties of the operator 
d xxxx . On the contrary, for general values of k, all critical values "moves" when k increases, 
as well as in [7]. Hence, it may happen that two different bifurcation values overlap for 
special values of k, in which case they are referred as resonant values. 
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Figure 1. The bifurcation picture for k — 1. 



Assuming that k = 1, for istance, Fig. 1 shows the bifurcation picture of solutions in 
dependence on the applied axial load p = bn 2 . In particular, u = and 

u f(x) = Af sin(Tra), A± = -y/2(b-l) , Af = ^2(b - 2) . 

4.1. The Lyapunov Functional and preliminary estimates. We begin to prove the 
existence of a Lyapunov functional for S(t), that is, a function C G C(7-L,M.) satisfying 
the following conditions: 

(i) C(z) — > oo if and only if \\z\\ n — > oo; 

(ii) £(S(t)z) is nonincreasing for any z eH; 

(iii) C(S(t)z) = C(z) for all t > implies that z G S. 

Proposition 4.3. // / is time-independent, the functional C defined in (3.2) is a Lya- 
punov functional for S(t) . 

Proof. Assertion (i) holds by the continuity of C and by means of the estimates 

~E(z)-c<£(z)<^E(z)+c. 
Using (3.3), we obtain quite directly 

(4.2) ^-aS(t)z) = -2\\d t u(t)\\ 2 <0, 

at 

which proves the decreasing monotonicity of £ along the trajectories departing from z. 
Finally, if C(S(t)z) is constant in time, we have that d t u = for all t, which implies that 
u(t) is constant. Hence, z = S(t)z = (u , 0) for all t, that is, z G S. □ 

The existence of a Lyapunov functional ensures that E(t) is bounded. In particular, 
bounded sets have bounded orbits. 

Notation. Till the end of the paper, Q : IRq — > IR + will denote a generic increasing 
monotone function depending explicity only on R and implicity on the structural constants 
of the problem. The actual expression of Q may change, even within the same line of a 
given equation. 



LONGTERM DYNAMICS OF THE EXTENSIBLE SUSPENSION BRIDGE 

Lemma 4.4. Given f G H , for allt > and initial data z e H with \\z\\ n < R, 
(4.3) £{t) < Q{R). 

Proof. Inequality (4.2) ensures that 

C{t) = C(S(t)z) < C(z) < Q(R) , Vt > 0. 

Moreover, taking into account that 



K0lr<^IK0ll2<^(*), 



we obtain the estimate 



C(t) > £(t) - 2 (/ , «(*)> > £(t) - 2 \\f\t 2 - \ \\u{t)f 2 > l -E{t) - 1 11 "' a 



Finally, we have 



£{t) < 2C(t) + 1 H/ll 2 < 2Q(i2) + 1 H/ll 2 = Q(i2). 



Lemma 4.5. Let p < v^i an d Fp(u) = Au — pA^u . Then 

(T p {u) , u) > C(p) \\u\\l, 

where 

'l, p<0 

(4.4) C(p) = " 

Proof. Because of the identity 



(i-^r), o<p <v ^. 



(J>(w), u) = \\u\\ 2 2 - P\\u\\l, 
the thesis is trivial when p < 0. On the other hand, when < p < y/Xi we have 

|2 



(J>(w) , u) = ||u||2 - p ||m|| 2 > ^1 - -^=^ 



l M ll2- 



We are now in a position to prove the following 
Theorem 4.6. When f = 0, the solutions to (1.3) -(1.5) decay exponentially, i.t 

£{t) < c £(0)e- ct 
with Co and c suitable positive constants, provided that p < y/X\. 
Proof. Let $ be the functional 

^( z )=S(z)+eT(z)- 1 -p 2 , 

where the constant 

(4.5) e = mm{\ 1 C(p),l} 
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is positive provided that p < V%~- In view of applying Lemma 4.5, we remark that 

$ = (F p {u) , u) + \\d t u\\ 2 + \ \\u\\\ + e (u, d t u) + k 2 \\u + \\ 2 . 
The first step is to prove the equivalence between £ and $, that is 

(4.6) ^£<*<Q{M\h)£- 

By virtue of (2.1), (4.5) and Lemma 4.5 the lower bound is provided by 

On the other hand, by applying Young inequality and using (2.1), we can write the 
following chain of inequalities which gives the upper bound of $. 

* < (c(p) + £ + 5 i-)|| lI |l5+(l + |)||9, ll || 2 -p|| ll |i; + i|| lI |lt< 
In particular, from (4.3) and (4.4) we find 



k 2 1 e 2 Q(\\z 



The last step is to prove the exponential decay of $. To this aim, we obtain the identity 

^ $ + £ $ + 2 (l - £ ) \\d t u\\ 2 + £ - \\u\\\ + e (1 - e) (d t u , u) = 0, 
where e is given by (4.5). Exploiting the Young inequality and (4.6), we have 

-$ + £$+(!-£) \\d t u\\ < ' \\u\\l < -^-^ $ , 



from which it follows 



d _, e(l + e) _ ^ n 

— $ H -$ < 0. 

dt 2 ~ 



Letting c = e (1 + e)/2, by virtue of Lemma 3.2 (with M — 0) and (4.6) we have 

2^ <*(*)<* (0) e' c * < £ (0) e~ ct . 

The thesis follows by putting c = 2 Ai<5(||^||«)/£. □ 

The existence of a Lyapunov functional, along with the fact that S is a bounded set, al- 
low us prove the existence of the attractor by showing a suitable (exponential) asymptotic 
compactness property of the semigroup, which will be obtained exploiting a particular de- 
composition of S(t) devised in [18] and following a general result (see [10], Lemma 4.3), 
tailored to our particular case. 
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4.2. The Semigroup Decomposition. By the interpolation inequality \\u\\l < ||m||||m||2 
and (2.1) it is clear that 

/ a rj\ 1 II i|2 ^ II ||2 II ||2 , || ||2 ^ || ||2 

(4.7) - ||tt|| 2 < 1 1 1 1 2 — + < m||u|| 2 , 

provided that a > is large enough and for some m = m(p, a) > 1. 

Again, i? > is fixed and ||z||% < -R. Choosing a > such that (4.7) holds, according 
to the scheme first proposed in [18], we decompose the solution S(t)z into the sum 

S(t)z = L(t)z + K(t)z, 

where 

L(t)z = (y(t), d t v{t)) and K(t)z = (w(t), d t w(t)) 
solve the systems 

d tt v + Av + <9 t -i; — [p — \\u\\\)A 1 l 2 v + av = 0, 
(v(0),d t v(0))=z, 



(4.8) 
and 
(4.9) 



d tt w + Aw + <9 t u> — (p — llull^A 1 / 2 ^ — av + k 2 u + = f, 
{w(0),d t w(0)) = 0. 



The next three lemmas show the asymptotic smoothing property of S(t), for initial data 
bounded by R. We begin to prove the exponential decay of L(t)z. Then, we prove the 
asymptotic smoothing property of K(t). 

Lemma 4.7. There is u = oo(R) > such that 

\\L(t)z\\ n < Ce-"K 

Proof. After denoting 

£o{t) = £ (L(t)z) = \\L(t)z\\ 2 n = \\v(t)\\l + \\d t v(t)\\ 2 , 
we set $o(*) = ®o(L(t)z,u(t)), where u(t) is the first component of S(t)z and 

«(*)) = S>(£(f)*) - pII II? + (« + \) Ht) ir + IKOII?IK*)ll? + <^(*),«(*)>- 

In light of Lemma 4.4 and inequalities (4.7), we have the bounds 
(4.10) l -E Q < $ < Q(R)£q. 

Now, we compute the time-derivative of $ along the solutions to system (4.8) and we 
obtain 

+ $ = 2(d t u,A 1 / 2 u)\\v\\ 2 1 < Q(R)\\d t u\\$ . 

The exponential decay of $ is entailed by exploiting the following Lemma 4.8 and then 
applying Lemma 6.2 of [18]. From (4.10) the desired decay of £ follows. □ 
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Lemma 4.8. For any e > 

t 

J \\d t u(s)\\ds<s(t-r) + £ l + ^p-, 

T 

for every t > r > 0. 

Proof. After integrating (4.2) over (r, t) and taking (4.3) into account, we obtain 

t 

-£{S{t)z) - 2^p- < + 2 [ \\d t u (s) f rfs = £(S(r)z) < £(*). 

2 Ai y 



It follows 



||$« (s)|| 2 rfs<g(i?), 



which, tanks to the Holder inequality, yields 



t 



J \\d t u(s)\\ds<eVt^+^p-<e(t-r) + ^ ' Q{R) 



e 



for any e > 0. □ 

The next result provides the boundedness of K(t)z in a more regular space. 
Lemma 4.9. (see [18], Lemma 6.3) The estimate 

\\K(t)z\\ n2 < Q(R) 

holds for every t > 0. 

Proof. As well as in [18], we use here the interpolation inequality 

HHI3 — IIHhllHU- 

Jointly with ||w|| 2 < Q(R) (which follows by comparison from (4.3) and Lemma 4.7), this 
entails 

(4.11) pWwgK^ + QiR), 
where 

S l (t)=S 1 (K(t)z) = \\K(t)z\\ 2 H2 = \\w(t)\\l + \\d t w(t)\\l 

Letting 

$i = £i + (\\u\H - p)\\w\\t + (d t w, Aw) - 2(f, Aw) + 2k 2 (u + , Aw), 
we have the bounds 

(4.12) ^£i-Q(R) <$i<Q(R)£i + Q(R). 
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Taking the time- derivative of $ l5 we find 

+ $1 = 2(d t u, A l ' 2 u) \\w\\l + 2a(A l ' 2 v, A 1 ' 2 d t w) + 



+ 



a(A 1/2 v, A 1/2 w) - (f, Aw)] - k 2 (Aw, u + ) + 2k 2 (Aw, d t u + ) 



Using (4.3) and (4.11), we control the rhs by 

+ Q{R)^S l + Q(R) < \ex + Q(R) < ^ + Q(R), 

and we obtain 

Since $i(0) = 0, the standard Gronwall lemma yields the boundedness of $i. Then, by 
virtue of (4.12), we obtain the desired estimate for E\. □ 

By collecting previous results, Lemma 4.3 in [10] can be applied to obtain the existence 
of the attractor A and its regularity. Within our hypotheses and by virtue of the decompo- 
sition (4.8)-(4.9), it is also possible to prove the existence of regular exponential attractors 
for S(t) with finite fractal dimension in T-L. This can be done by a procedure very close to 
that followed in [18]. Since the global attractor is the minimal closed attracting set, we 
can conclude that the fractal dimension of A in Ti is finite as well. 
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